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A combination of neutron total scattering measurement and reverse Monte Carlo (RMC) re- 
finement is applied to the study of apical Cu-O bond distortions in the high-Tc superconductor 
YBagCuaOe.ga- We show that the average structure is not consistent with a split-site model for the 
corresponding Cu and O positions, but that the local structure nevertheless reveals the existence of 
two separate apical Cu-O bond lengths. Using G{r) data obtained from a variety of Qmax values 
we show that this result is independent of the data treatment methodology. We also find that the 
resulting 'short' and 'long' Cu-O bond lengths agree well with the results of previous EXAFS stud- 
ies. The existence of bimodal apical Cu-O bond distributions in the context of a single-site average 
structure model is interpreted in terms of correlated displacements of the Cu and O atoms. We find 
evidence also for the clustering of short apical Cu-O bonds within our RMC configurations. 



I. INTRODUCTION 



The mechanistic role of structural inhomogeneities 
has remained a contentious aspeclPI^ of the sci- 
ence of high-temperature superconductors almost since 
the phenomenon was first discovered — be it in the 
form of electron-phonon coupling,'^ charge localisation 
('stripes'^P or bipolaron formation and condensation.'^ 
Local-structure probes, such as X-ray absorption spec- 
troscopy and total scattering (or pair distribution func- 
tion, PDF) measurements, might be expected to be par- 
ticularly sensitive to the extent and nature of lattice 
distortions; yet the results of such experiments have of- 
ten appeared contradictory and remain controversial. As 
a consequence there exists no generally-accepted micro- 
scopic description of local structure distortions in canoni- 
cal high-Tc materials such as the YBa2Cu307_5 (YBGO) 
family. 

In the particular case of YBGO, it is the existence 
of a double-well potential for the apical Gu2-04 bond 
that has proven particularly controversial, primarily be- 
cause of the associated implications for charge localisa- 
tion and electron-lattice coupling [Fig. [l]ja)].^ An early 
XANES study was the first to suggest the existence of 
such a double-well,-^ and a number of EXAFS inves- 
tigations in the intervening years have essentially sup- 
ported th is sam e conclusion.lSHSl in contrast, both sin- 
gle crystaP^^ and powdei^^^ltlSl diffraction experiments 
found no evidence for splitting of either the Gu2 or 04 
crystallographic sites. Likewise neutron PDF studies 
ruled out a bimodal distribution of 04 positions but 
yielded improved fits for models tha t allo wed a variety 
of different splittings of the Gu2 site.'^SHH And whereas 
frozen-phonon ab initio calculations failed to reproduce 
a double-well potential for the Cu2-04 bond,!^ inelas- 
tic neutron scattering measurements linked an increased 
sensitivity of the dynamic scattering function S{Q, E) 
near Tc to structural instabilities involving shortening of 




FIG. 1: Representations of the average structure of 
YBa2Cu3 06.93 determined by (a) Rietveld refinement and (b) 
RMC refinement of neutron total scattering data. The crucial 
Cu2-04 bond is indicated by the arrow. Thermal ellipsoids 
in (a) are shown at 80% probability level; the representation 
in (b) is obtained by projecting the RMC configuration onto 
a single unit cell. 



the same bond.'^^l 

In this paper we attempt to rationalise some of these 
apparently contradictory results by carrying out a re- 
verse Monte Carlo (RMG) refinement of newly-collected 
neutron total scattering data for a carefully synthesised 
YBGO sample of composition YBa2Gu306.93. This is 
very much a PDF study, but there are two key differ- 
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ences between our approach and those of Refs. [2D] and 
[2T] First, where these previous studies have focussed on 
interpreting structural features of the PDF in terms of in- 
creasingly complex single-particle correlations — through 
symmetry lowering of the unit cell or the introduction 
of split sites — RMC refinement allows the PDF to be in- 
terpreted in terms of correlated displacements of atoms. 
Second, our RMC refinements include simultaneous fits 
to both PDF and Bragg intensity data, a process that en- 
sures the resultant atomistic configurations are automati- 
cally consistent with traditional average-structure refine- 
ments. The key result of our study is that the Cu2-04 
bond distribution in our RMC configurations is indeed 
bimodal, but that correlations between the Cu2 and 04 
displacements result in unimodal distributions for both 
Cu2 and 04 sites. In this way, we are able to rationalise 
the apparent discrepancy amongst EXAFS, PDF and av- 
erage structure crystallographic studies. 



II. MATERIALS AND METHODS 

A. Synthesis and characterisation 

A sample of YBa2Cu307_5 (2.5 g) was prepared via the 
citrate gel route described previously.!^ A suitable stoi- 
chiometric mixture of Y2O3 (99.999%, dried at 900 °C), 
BaCOa (99.997%) and CuO (99.9999%) was dissolved 
in a minimum quantity of 6 M nitric acid. Four mole 
equivalents of citric acid and 4 ml of ethylene glycol were 
then added and the solution heated with constant stir- 
ring. The gel formed was subsequently ground into a 
fine powder, placed in a crucible lined with Ni foil, and 
heated in air to 600 °C at a rate of 1°C min^^. The 
resulting powder was re-ground and pressed into 13 mm 
pellets under a force of 5t. All subsequent heating was 
performed with the sample contained within an alumina 
boat lined with Ni foil. The sample was fired at 920 °C 
under flowing O2 and then cooled to room temperature 
at a rate of l°Cmin~^. This process was iterated three 
times. After the final cooling cycle, the pellets were re- 
ground and the resulting powder heated at 600 °C for 12 h 
under flowing O2 to give a sample with O substoichiom- 
etry 6 = 0.07(3) as determined by iodometric analysis. 
Phase purity was confirmed by X-ray powder diffraction; 
the refined lattice parameters were found to be in good 
agreement with previously published values. A mea- 
surement of the superconducting transition temperature 
by SQUID magnetometry gave = 92(1) K. 



B. Neutron total scattering 

Neutron total scattering data were collected using 
the GEM instrument at ISIS.'2£H2l A sample of poly- 
crystalline YBa2Cu3 06.93 (1-8 g), prepared as described 
above, was placed within a cylindrical thin-walled vana- 
dium can of 3 mm diameter and 5.8 cm height, which was 



in turn loaded inside a closed cycle helium refrigerator. 
The sample was cooled to 50 K and total scattering data 
then collected over a large range of scattering vectors of 
magnitudes 0.7 < Q < 50 A~^, corresponding to a real- 
space resolution of order Ar ~ 3.791/Qmax — 0.08 A. 

Following their collection, the total scattering data 
were corrected using standard methods, taking into 
account the effects of background scattering, absorp- 
tion, multiple scattering within the sample, beam inten- 
sity variations, and the Placzek inelasticity correction.l^Sl 
These corrected data were t hen co nverted to experimen- 
tal F{Q) and G'(r) functionsPSSI 



F{Q) = po I 4^r2G(r)=^dr 



(1) 



G(r) = c,c,&,6,fe(r)-l], (2) 



where 



Airr'^drpj 



(3) 



riijir) the number of pairs of atoms of type i and j sep- 
arated by distance r, po is the number density, Ci the 
concentration of each species i and bi the corresponding 
neutron scattering length. As part of the data normali- 
sation process it is usual practice to determine a useable 
value Qmax that represents the best possible compromise 
between optimising Ar for PDF refinement {i.e. mak- 
ing Qmax as large as possible) and avoiding the unnec- 
essary inclusion of high-frequency noise that can plague 
high-Q data. In this instance, we found that a value 
Qmax — 40 A^^ produced the most refiable G(r) func- 
tion. 

The Bragg profile functions for each data set were ex- 
tracted from the scattering data collected by the detector 
banks centred on scattering angles 29 ~ 54.46°, 63.62° 
and 91.37°. The experimental Bragg diffraction profiles 
were fitted with the GSAS Rietveld refinement progranPS 
using the published structural modelp^l the fitting pro- 
cess is discussed in greater detail in the Results section 
below. 



C. Reverse Monte Carlo refinement 

The reverse Monte Carlo refinement method as applied 
to crystalline materials, together with its implementa- 
tion in the progra m rm cprofile have been described 
in detail elsewhere .^^SEH The basic refinement objective 
is to produce large atomistic configurations that can ac- 
count simultaneously for the experimental F{Q), G{r) 
and Bragg profile I{t) functions. This is achieved by 
accepting or rejecting random atomic moves subject to 
the metropolis Monte Carlo algorithm, where in this case 
the Monte Carlo acceptance criterion is determined by 
the quality of the fits to data. The refinement process 
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TABLE I: 'Distance window' parameters dmin, rfmax used for 
all RMC refinements in this study. 



Atom pair 


^^min (A) 


dmax. (A) 


Y-02,3 


2.f0 


2.56 


Ba-O 


2.50 


3.50 


Cul-Of 


1.73 


2.10 


Cul-04 


1.73 


2.50 


Cu2-02 


1.73 


2.10 


Cu2-03 


1.73 


2.10 


Cu2-04 


1.73 


2.50 



is continued until no further improvements in the fits to 
the data are observed. By virtue of the particular impor- 
tance in the present study of fitting accurately the lowest- 
r region of the PDF — which includes the all-important 
distribution of Cu2~04 separations — we applied a larger 
weighting to the region < r < 7A than to the remain- 
der of the PDF. 

Our starting configurations for the RMC process were 
based on a 24 x 24 x 8 supercell of the crystallographic 
unit cell shown in Fig.jlja). Each configuration contained 
59 904 atoms and extended approximately 90 A in each 
direction. By virtue of the small experimental value of S 
for our sample, and in order to avoid introducing model 
bias, we did not incorporate any explicit consideration 
of O vacancies within our model. In addition to the re- 
quired fits to data, the only constraints placed on the 
atomic coordinates were a set of 'distance window' con- 
straints which act to maintain an appropriate framework 
connectivity throughout the refinement process.''"^ "^^ In 
each case the values used [see Table |l] were based on the 
extrema of the corresponding peaks in the experimental 
PDF. 



III. RESULTS 



TABLE II: Crystallographic parameters, atomic coordinates 
and isotropic equivalent displacement parameters determined 
using Rietveld refinement of neutron scattering data for 
YBaaCusOe.gs. 



Crystal system Orthorhombic 

Space group Pmmm 

a (A) 3.81412(5) 

b (A) 3.87694(6) 

c (A) 11.63970(22) 

V (A^) 172.1170(30) 
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2 


f/iso (A2) 


Y 


0.5 


0.5 


0.5 


0.0031(4) 


Ba 


0.5 


0.5 


0.18377(12) 


0.00092(26) 


Cul 











0.0038(4) 


Cu2 








0.35535(9) 


0.0019(8) 


01 





0.5 





0.0115(7) 


02 


0.5 





0.37837(13) 


0.0045(4) 


03 





0.5 


0.37727(15) 


0.0020(4) 


04 








0.15879(11) 


0.0047(4) 



value was fixed at 0.001 A^. We remark also that dis- 
placements of the so-called 'chain' oxygen atoms (01) are 
clearly very anisotr opic. Th is behaviour, which has been 
reported elsewhere jlJ^^I^IESI jg thought to be related to 
low-energy correlated displacements of the Cul-centred 
[CUO4] square units. 

Since our study concerns primarily the distribution of 
Cu2-04 bonds, we investigated the validity of structural 
models which incorporated splitting of one or other or 
both of the two atomic sites. As in previous Rietveld 
studies,'^!^^^ we found that none of these three cases 
yielded a stable refinement. 



A. Average structure determination 

Our first step was to use the GSAS refinement packag^^ 
to verify that the Bragg contribution to our total scatter- 
ing data could be interpreted in terms of sensible lattice 
parameters, atomic coordinates and anisotropic displace- 
ment parameters. The values obtained, which are sum- 
marised in Table |TIj are in good agreement with previous 
crystall ographic studies of nearly-stoichiometric YBCO 
samples .ISHSMIrpj^g corresponding fit to data is shown in 
Fig. [2ja) and a representation of the structural model 
obtained is that shown in Fig. [ija). It was possible 
to achieve stable refinement of anisotropic displacement 
parameters for all atoms other than Ba, for which un- 
constrained refinement consistently yielded non-positive 
definite values for the U22 parameter; consequently this 



B. Local structure determination 

RMC refinement yielded equally satisfactory fits to the 
neutron scattering data — both in terms of the reciprocal- 
space Bragg intensity function I{t) and the real-space 
G{r) transform. Importantly, the quality of the Bragg 
profile fit is comparable to that obtained from Rietveld 
refinement, implying that the RMC configurations are as 
consistent with the Bragg component to the total scat- 
tering function as the average-structure model given in 
Table HH 

By collapsing the atomic coordinates of an entire RMC 
configuration onto a single unit cell, it is possible to vi- 
sualise the RMC-refined average structure [Fig. [ijb)]. 
As found in the Rietveld model, the 01 distribution is 
anisotropic (and indeed there is some sign here of a split 
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FIG. 2: The fit to neutron powder difTraction data obtained 
by (a) Rietveld and (b) RMC refinement, (c) The correspond- 
ing G(r) fit obtained using RMC. In all cases, data are shown 
as filled circles, fits as red lines and the difference (data— fit) 
as blue lines, sometimes shifted vertically for clarity. A re- 
fined background function has been subtracted from the data 
and fit shown in (a). 



site, although we do not pursue this possibility further 
here). Encouragingly, the Ba site is now represented 
by a physically sensible distribution function; because 
RMC configurations represent a physical box of atoms, 
there is never an issue of obtaining non-positive-definite 
anisotropic displacement parameters. 

Because the RMC-refined distributions of atomic posi- 
tions are not constrained to assume ellipsoidal forms, we 
are able to check in a straightforward and unbiased man- 
ner whether or not there is any real evidence for splitting 
on either the Cu2 or 04 sites. In both instances a projec- 
tion of the corresponding atomic coordinates onto the z 
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FIG. 3: (a) Cu2-04 partial pair distribution functions (solid 
points) determined by RMC refinements against G{r) func- 
tions generated for maximum scattering vectors of (left) 
Qmax = 40 A"^ and (right) Qmax = 30A"\ These two dis- 
tributions correspond respectively to the most and the least 
convincing bimodal distributions obtained for the 11 differ- 
ent Qmax values investigated. Calculated fits using a pair of 
Gaussian curves are shown as a solid line, with the two indi- 
vidual Gaussian contributions shown as red and blue dashed 
lines, (b) Qmax-dependence of the 'short' and 'long' Cu2-04 
bond lengths (red and blue circles, respectively), determined 
as the midpoints of the corresponding Gaussian fits. Error 
weighted averages are plotted as vertical red and blue dashed 
lines. The range of values obtained using EXAFS in Refs. ITOl 
and [12] are shown as shaded regions, with the Rietveld Cu2- 
04 bond length value shown as a bold black vertical line. 



axis gave distributions that could be well fitted using sin- 
gle Gaussian functions. Hence we can say that our RMC 
refinements are consistent with unimodal single-particle 
correlation functions for both atom sites. In particular, 
there can be no need to invoke split Cu2/04 atom sites in 
order to fit the neutron scattering data since RMC would 
have been free to do so should the data have demanded. 
With the RMC description of average structure in 
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YBa2Cu306.93 established, we turned our attention to 
the pair correlations at the heart of the local structure 
controversy. It is straightforward to calculate from an 
RMC configuration the distribution functions for specific 
bonding interactions. Carrying out such a calculation 
for the crucial apical Cu2--04 bonds yields a distribu- 
tion that is quite obviously bimodal [Fig. |3ja) , left-hand 
panel]. The distribution is will fitted using a sum of 
two Gaussian components, ascribable to contributions 
from 'short' and 'long' Cu2-04 bonds. The midpoints 
of these Gaussian fits suggested a pair of Cu2-04 bond 
lengths that were surprisingly similar to the values ob- 
tained in the EXAFS study of Ref. IH rcu2-04 = 2.167 
and 2.307 A (RMC) vs rcu2-04 = 2.220 and 2.337 A (EX- 
AFS). 

We were interested to establish the origin of this fea- 
ture in the experimental G(r) function in order to de- 
termine the extent to which the close agreement with 
EXAFS results is fortuitous. There is a clear peak in the 
G(r) function responsible for splitting the Cu2-04 dis- 
tribution [Fig. |4|a)]. The only non-zero partial pair dis- 
tribution function in this particular region is 5cu2,04('')- 
All other Cu-0 bonds are much shorter and contribute 
to the first significant peak seen at r ~ 1.9 A. There is 
some overlap at higher values of r between the Cu2-04 
and Y-02/Y-03 partial pair distribution functions, but 
the values of gY,02ir) and gY,03('') are essentially zero 
for r < 2.20 A, despite the RMC constraints allowing Y- 
O distances to be as small as 2.10 A. Consequently the 
r ~ 2.15 A feature in the G(r) function — if real — can be 
ascribed only to a separate distribution of short Cu2-04 
bonds, and hence is an experimental indication of the 
existence of a bimodal Cu2-04 bond distribution. 

However, the key question here is whether or not the 
feature is indeed real — our concern being that its magni- 
tude is comparable to the size of the Fourier truncation 
ripples evident at the lowest values of r. Consequently, 
we proceeded to check whether or not the existence of a 
bimodal distribution was affected by the value of Qmax 
used to generate the G(r) function. Changing Qmax 
varies strongly the location and periodicity of Fourier rip- 
ples but ought not greatly affect the positions and intensi- 
ties of any real features of the data (though one must keep 
in mind that reducing Qmax also reduces the real-space 
resolution Ar). To this end we calculated G(r) data for 
11 equally-spaced values 30 < Qmax < 50 A"^ [Fig.[4];b)] 
and used each of these modified data sets to drive an 
independent RMC refinement. On convergence, refine- 
ment was continued for sufficiently many RMC steps to 
obtain a minimum of seven independent configurations 
for each run (independent configurations being separated 
by 140 000 ~ 59 904 In 10 accepted movePl). 

For each of the 11 sets of RMC configurations pro- 
duced in these refinements, the corresponding Cu2-04 
bond-length distributions were calculated and fitted as 
described above. It was found that two Gaussian func- 
tions were always required for these fits and that the 
corresponding Gaussian midpoints did not vary particu- 
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FIG. 4: (a) Experimental G(r) function (filled circles) and 
corresponding RMC fit (solid black line). The relative contri- 
butions to the RMC fit arising from the various partial pair 
distribution functions are included, with those corresponding 
to the Cu-O, Y-0 and Ba-O pairs highlighted in red, green 
and blue, respectively. The experimental feature responsible 
for spHtting the Cu2-04 distribution is highlighted with an 
asterisk, (b) Experimental G(r) functions generated using a 
range of Qmax values 30 < Qmax < 50 A~^; the bold curve is 
that shown in (a). 



larly strongly with Qmax [Fig. |3][c)]. The values of the 
short and long Cu2-04 bond lengths we obtain as aver- 
ages over the various configurations are compared with 
the EXAFS results of Ref. [H] in Table [ml To give an idea 
of the level of confidence in these values, we have shown 
examples of the distributions corresponding to both the 
most and the least convincing bimodal fits [Fig. [3]ja)]. 
There is an additional consistency in the Gaussian widths 
and peak areas: in both cases the values obtained across 
all 11 data sets vary by less than 10%. Taking an error- 
weighted average across all 11 sets of configurations, we 
obtain a fraction of short Cu2-04 bonds of 12.3±2.8%. 
For completeness, we note that it is not at all unexpected 
that the bond lengths extracted for the lowest and high- 
est Qmax values are the least reliable: in the former case, 
the real-space resolution Ar is heavily reduced, and in 
the latter case the incorporation of high-Q noise begins 
to affect noticeably the smoothness of the G(r) function. 
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TABLE III: Free atom position parameters and Cu2-04 bond lengths determined using various experimental techniques and 
refinement strategies. 
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ZCu2 


0.35535(9) 


0.35464(8) 


0.3570(49) 


0.3548(2) 


0.35580(7) 




Z02 


0.37837(13) 


0.37808(7) 


0.3778(29) 


0.3780(3) 


0.37701(11) 




Z03 


0.37727(15) 


0.37808(7) 


0.378(2) 


0.3782(3) 


0.37805(8) 




204 


0.15879(11) 


0.15919(10) 


0.15998(95) 


0.1599(2) 


0.15951(9) 




(d(Cu2-04)) 


2.2879(17) 


2.2860(15) 


2.268(14) 


2.280(3) 


2.2849(13) 




d(Cu2-04),hort 






2.20(10) 


2.19(3) 


2.197(10) 


2.220(5) 


d(Cu2-04)long 






2.29(10) 


2.37(3) 


2.306(3) 


2.337(5) 



C. PDFgui refinements 



In order to ensure consistency with previous PDF stud- 
ies, we performed a final set of PDF refinements using the 
PDFGUI software package.^^ Sometimes termed a 'real- 
space Rietveld' approach, PDFGUI uses the experimental 
G(r) function to refine atomic coordinates and displace- 
ment parameters within a single unit cell. Our aim here 
was to reproduce the results of Ref. i21j namely, that a 
split site could be refined for the Cu2 atom but not for the 
apical 04 atom, and that the split Cu2 atom sites corre- 
sponded to a difference A(i(Cu2-04) between 'short' and 
'long' Cu2-04 bond lengths of 0.18(6) A. 

Our PDFGUI refinements yielded the fit shown in Fig. [5] 
where we have used a Qmax value of 25 in order to 
replicate the conditions of Ref. ,21j It was noted in this 
earlier study that the range of r values included in the 
fitting process had a significant effect on the value of 
Ac?(Cu2-04) obtained; this r-dependency in PDFGUI re- 
finements is symptomatic of the existence of short-range 
correlations that are not well described by the long-range 
structural periodicity (see e.g. Ref. 36]. Consequently, 
we determined values of Ad(Cu2-04) for the two values 
employed in Ref. (rmax = 5 and 15 A) and also for 
a range of G(r) functions generated using increasingly 
large Qmax values; our results are given in Tables |III| and 

CYl 

While we found that it was possible to refine split 
sites for the Cu2 atom, the difference between 'long' and 
'short' Cu2-04 bonds was slightly smaller than that re- 
ported in Ref. |21j Like the earlier study, we also found 
that split-site refinements gave improved fits when only 
the lowest-r region of the PDF was used, suggesting that 
the effects of site splitting only impact noticeably on the 
PDF only over short distances. We note that, for some 
of the various combinations of Qmax and rmax, it was ac- 
tually possible to refine a split 04 site; however the cor- 
responding distance between sites was never larger than 
0.05 A. 




r(A) 

FIG. 5: A typical fit to G(r) data obtained using PDFGUI, 
in this instance for the experimental G(r) function obtained 

using Qmax = 25 A"^ 



TABLE IV: Differences Ad(Cu2-04) between short and long 
Cu2-04 bond lengths determined using r-dependent PDFGUI 
refinement of the experimental G(r) functions generated for 
various Qmax values. 


(A-) 


(A) 


Ad(Cu2-04) 

(A) 


Rws 


Rwu 


25 


5 


0.08(7) 


0.0908 


0.0935 




15 


0.03(5) 


0.0883 


0.0901 


30 


5 


0.089(26) 


0.1224 


0.1406 




15 


0.03(6) 


0.0941 


0.0951 


40 


5 


0.099(28) 


0.1207 


0.1267 




15 


0.065(26) 


0.1261 


0.1276 


50 


5 


0.093(22) 


0.1863 


0.1888 




15 


0.008(4) 


0.2146 


0.2157 



IV. DISCUSSION 

Perhaps the key result of our RMC study is to demon- 
strate that a single atomistic configuration can account 
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at once for a number of the seemingly-disparate experi- 
mental results reported previously forYBa2Cu307_5: 

1. The finding h'om both single crystal and powder 
diffraction studies that the scattering distributions 
at the Cu2 and 04 sites of the average structure 
are unimodal. 

2. That two Cu2-04 bond lengths are required to ob- 
tain satisfactory fits to Cu K-edge EXAFS spectra. 

3. For PDFGUI refinements over small regions in real 
space (r^ax — SA"^), neutron PDF data are best 
fitted in terms of a structural model with two Cu2 
sites (and hence two Cu2-04 distances); however 
as Tmax is increased the refined splitting between 
Cu2 sites vanishes. 

The clear implication is that these results are not con- 
tradictory, but merely reflect the importance of taking 
into account length scale when interpreting structural 
features of YBa2Cu307_5. 

Consistency between a bimodal Cu2-04 bond length 
distribution on the one hand, and single Cu2/04 sites 
in the average structure on the other hand, demands 
that the Cu2 and 04 displacements are correlated during 
short- and long-bond formation. Specifically, shortening 
of the Cu2-04 bond must involve cooperative displace- 
ment of the Cu2 atom in the —z direction [with reference 
to the coordinates and axes of Fig. [ija)] and of the 04 
atom in the +z direction; lengthening of the same bond 
must involve correlated displacement in the opposite di- 
rection. A split-site model for Cu2 and 04 positions 
does not capture this situation properly in either average- 
structure or PDF refinements. In the former case the 
model fails because displacements away from the average 
Cu2/04 sites will be indistinguishable from thermal mo- 
tion if correlation is not taken into accountp^ in the lat- 
ter case — and again because correlation is ignored — the 
model actually predicts a trimodal Cu2-04 distribution, 
which is clearly inconsistent with the bimodal distribu- 
tion observed experimentally. 

For the low-r„iax PDFGUI refinements, splitting of the 
Cu2 site allows the bimodal Cu2-04 bond distribution 
to be fitted, but does so with a model that is incompat- 
ible with the average structure; hence if the same model 
is fitted over larger regions in real-space then either the 
fit is degraded or, if the magnitude of the Cu2 splitting 
is allowed to refine, a single-site model is recovered. The 
key advantage of supercell refinement techniques such as 
RMC is that the large box size allows sufficient sampling 
to capture both local and periodic structural features 
within a single atomic-scale model. 

Certainly our RMC refinements suggest that in 
YBa2Cu3 06.93 there are two types of apical Cu2-04 
bond and, as already discussed, the difference in bond 
length that we measure (0.11 A) is remarkably consis- 
tent with the results of previous EXAFS studies. There 
will of course be a temptation to associate with the two 
Cu2-04 bond types a different valence state for the Cu2 
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FIG. 6: The short-bond/short-bond correlation function cal- 
culated for one of the RMC configurations described in the 
text. For ease of interpretation, the horizontal axis corre- 
sponds to interatomic distances scaled relative to the average 
in-plane lattice parameter. The large positive peak at r/ucff 
indicates that short apical Cu2-04 bonds are likely to form 
neighbouring pairs. 



atom; however we note only that the neutron PDF mea- 
surements we have performed here involve no scattering 
contribution from the electronic structure of the material 
and hence we do not comment further on this aspect. 

We might, however, expect our RMC configurations to 
be sensitive to correlations in the spatial arrangement of 
the different types of Cu2-04 bonds. We have been made 
aware of an earlier, unpublished, RMC study in which 
short-bond/short-bond correlations are discussed,!^ and 
we conclude here with a similar calculation for our own 
RMC configurations. In order to quantify the extent to 
which short Cu2-04 bonds are found to cluster, we define 
a correlation function 

i 3 

where 

[ +1 d(Cu2-04)i > 4 
Ui = { , (5) 

i -1 d(Cu2-04)i < 4 

the distance dc is the critical Cu2-04 bond length below 
which a bond is considered to be 'short', N is the num- 
ber of Cu2 centres in a configuration, and the sums in 
Eq. |4]are taken over all Cu2 centres i,j within the same 
[CUO2] layer and separated by the distance r. Because 
the number of short and long bonds differ, x(^) '^iU gen- 
erally be non-zero even for random distributions of short 
bonds; consequently we calculate also a statistical cor- 
relation function Xrand('') which is generated using the 
same set of Ui values as x(^)i but with these distributed 
randomly throughout the configuration. Then the differ- 
ence function 

Ax(f) = X{r) - Xrand(7-) (6) 
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measures the extent to which short Cu2-04 bonds are 
more (Ax(r) > 0) or less {Ax{r) < 0) Ukely to cluster at 
a distance r than in a random distribution. 

Calculation of Ax{r) for any one of our RMC config- 
urations yields a correlation function of the general form 
shown in Fig. [6j In all cases we find only one consistently 
well-defined peak: namely a positive value of A.x{r) at 
r = tteff (= Vab). This would suggest that, in general, 
short Cu2-04 bonds are more likely to cluster into pairs 
than to be distributed statistically throughout the config- 
uration. If accompanied by charge localisation, then such 
a result could be considered consistent with a bipolaron 
modellSl 

In conclusion, RMC refinements of neutron PDF data 
collected for a YBa2Cu306.93 sample yield atomistic con- 
figurations that are simultaneously consistent with all 
features of the short-range and long-range structural fea- 
tures captured by total scattering experiments and with 



the results of independent X-ray absorption spectroscopy 
experiments. We find no evidence to support the use of 
a split site model, but do find that the apical Cu2-04 
bonds are of two types: a short bond of 2.197(10) A and 
a long bond of 2.306(3) A. There appears to be some clus- 
tering of the short Cu2-04 bonds into nearest-neighbour 
pairs within [CUO2] layers. 
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